Lattice Induced Resonances in One Dimensional Bosonic Systems 
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We study the resonant effects produced when a Feshbach dimer crosses a scattering continuum band of atoms 
in an optical lattice. We numerically obtain the exact spectrum of two particles in a one-dimensional lattice 
and develop an effective atom-dimer Hamiltonian that accurately captures resonant effects. The lattice-induced 
resonances lead to the formation of bound states simultaneously above and below the scattering continuum 
and significantly modify the curvature of the dimer dispersion relation. The nature of the atom-dimer coupling 
depends strongly on the parity of the dimer state leading to a novel coupling in the case of negative parity dimers. 

PACS numbers: 



Feshbach resonances HI and optical lattices |2| are two 
powerful experimental tools used to drive ultracold atomic 
gases into the strongly correlated many-body regime. While 
recent experiments have realized such resonant lattice sys- 
tems 1 3 1, current understanding of them is limited, with tra- 
ditional single-band Hubbard model failing to capture the res- 
onant regime. However, a many-band description necessary 
near a resonance (where interaction is much larger than the 
band width and gap), when implemented directly may be un- 
necessarily complicated and impractical for analytical stud- 
ies. Thus, the development of better suited effective models 
is highly desired and important first steps have been taken in 
this direction HUS). 

At the two-body level, the lattice potential changes the scat- 
tering properties, shifting existing resonances and inducing 
new ones. This occur when dimers formed with atoms in 
excited bands cross and hybridize with the scattering contin- 
uum of the lowest band |6|. While resembling confinement- 
induced resonances |7|, lattice resonances exhibit richer phe- 
nomenology and are more challenging to describe due to the 
coupling between the relative and the center of mass degrees 
of freedom. Previous two-body lattice studies [8| have mainly 
focused on the lowest molecular state and the lowest atomic 
band with the virtual higher renormalizing the parameters of 
the single-band Hubbard model. 

In this Letter, we study the ID lattice-induced resonances. 
We derive an effective atom-dimer many-body model that cap- 
tures all the features of our exact numerical solution of the 
two-body problem summarized in Figs. [T] and |3] These in- 
clude: (i) attractive and repulsive resonance-induced bound 
states, that for a range of the center-of-mass momenta can si- 
multaneously appear above and below the two-particle contin- 
uum bands, (ii) ID lattice resonances (absent in ID lattice-free 
systems) at the center-of-mass momenta, where a bound state 
enters the two-particle continuum, (iii) strong interaction- 
dependence of the molecular bound-state dispersion, that can 
even be driven negative. Our two-body results are also appli- 
cable to two-species fermions in a singlet state. However, for 
many-system, in this Letter we focus on bosonic atoms. The 
effective Hamiltonian, inspired by Ref. [4], is a lattice pro- 
jection of a two-channel model |j9l- The novel feature of our 
model is a nearest neighbor atom (a^) - dimer (di) Feshbach 



coupling. 

Had ^ gY,dU^i^^+l + (-l)^"a«-i) + H.c. + (1) 

i 

that depends sensitively on the dimer's parity Pd ~ (sym- 
metric) Pd = I (antisymmetric), that is crucial to capturing 
the dispersion of lattice-induced resonances and bound states. 
We emphasize that the key atom-dimer coupling Had, Eq. ([T]) 
proposed here, describes two atoms forming a dimer via a pro- 
cess which depends on the dimer's momentum, thus is not re- 
lated to the usual two-channel model of Feshbach resonances. 
We derived Eq. ([TJ via a careful analysis of the dimer orbital 
structure, supported by exact diagonalization of a two-atom 
Hamiltonian [Eq. (|2]) below] in a periodic potential. 

Our starting point is the bosonic ID lattice Hamiltonian 
with contact interaction, 

H = ^ f - — + Vosm^{kiXa)j+^ gioSixa-xp), 

(2) 

where a is the lattice constant, ki = tt/o, and Vq is the lat- 
tice depth. The Hamiltonian (2) can be realized with current 
lattice experiments [10] with tight confinement a± <^ a\\ { 
a±,\\ are oscillator lengths lattice wells in the _L and || direc- 
tions) and negligible tunneling in the perpendicular direction 
( Jl ~ 0). The interactions are controlled by tuning the scat- 
tering length Qs close to a Feshbach or confinement-induced 
resonance |7|. In the regime |as| < a j^, the atoms remain well 
localized in the transversal ground state and the low-energy 
behavior is described by Eq. (|2]l with gio « Afi^as/imaY)- 
The typical lattice energy scale is Ej^ = h?kf /2m and the di- 
mensionless interaction strength is A = gioki/Er. In free ID 
space (Vb — 0), the two-body Hamiltonian supports a bound 
state at an arbitrarily weak interaction with a binding energy 
Eb - glom/iAh^). 

The two-boson Schrodinger equation is numeiically solved 
for a system of L sites and periodic boundary conditions us- 
ing a plane wave expansion ifTTj . The convergence has been 
extensively tested by changing L = 1, 21, and the basis 
dimension. Figures [T] and [3] summarizes the numerical results 
for different lattice depths and attractive interaction values. 
States outside the scattering continuum bands [colored regions 
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in Figs.[T](a) and|3) represent bound states, while the states in- 
side the continuum bands are either scattering states or dimer 
states which decoupled from the scattering continuum. 
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FIG. 1: (Color online) (a) Two-body spectrum for a lattice with 
Vo ~ 'iEr as a function of K. Colored regions correspond to the 
scattering continuum bands [lower region corresponds to Eq^ (green 
online), and the upper region corresponds to Ei^ (red online)]. The 
curves are the (0,0), (1,0) and (2,0) bound states at different in- 
teraction strengths A, On the lattice, the (2,0) dimer state is lower in 
energy than the (1,1) dimer state. The (2, 0) dimer state is lower in 
energy than the (1,1) dimer state, (b) and (c) Exact solutions (solid 
curves) and two-band Hubbard model predictions (dashed curves) of 
the energy (b) and effective tunneling J"-^-^ (c) for the lowest two 
bound states of a lattice of = lOEr- Shaded regions in (b) corre- 
spond to Eq'^q and E'^'^q and the shaded region in (c) to the interaction 
regime in which the excited dimer enters i?o,o- 



The single particle energies are grouped in to bands 
El^{k) with n > which in the tight binding regime 
take the form El^^ik) — ta.n — 2 Ja.„ cos(A;a), where 
^a,n — J \ wns{x)\'^ Ho{x)dx is the onsite energy and Ja,n — 
I ""^n i+i{x)Ho{x)wn,i{x)dx is the nearest neighbor tunnel- 
ing. Here, Hq is the noninteracting part of the H [2j and 
Wn.i{x) the wannier function of band n centered at site i lfT2l . 
The two-body solutions describe the scattering continuum 
bands [colored regions in Fig. [T] (a)] which are the sym- 
metrized product of the single-particle eigenstates with ener- 
gies El-„,{K,k) = El-{K/2 + fc) + £;^"(i^/2 - fc), and 
can be classified by the band label pair {n,m). For each 
scattering continuum, we can define the typical band gap as 
hujnm = min(|ea,„+ea,m-ea,fe-ea,/|) with {k, 1} ^ {n, m}. 

For weakly interacting atoms (|5i£)/a||| ^ hiOnm) in 
the tight-binding regime, a bound state is formed in the 
vicinity of each isolated (not overlapping) scattering con- 
tinuum {n,m) with parity Pd = (—1)"+™. A natural 
starting point to capture the molecular behavior in Fig. 1 is 
a tight binding model with each of the two particles with 
its own (possibly the same) hoping J„ and J^, and an 
on-site interaction — gio j \'Wm,i{x)\'^\wn,i{x)\'^dx 

(of order of gin/aw)- This model can be solved analytically 
giving a molecular spectrum E'j^^'^^^^ (K) = + fa,™ + 

sign{Unm) J U^m + ^Ja,n + 4Ja,m + ^Ja,nJa,m COs{Ka). 



Solid and dashed curves in Figs. [T| (b) and (c) 
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the lowest two scattering continua E'q'J, and E1°q. While the 
two-band description is accurate at weak interaction, it clearly 
breaks down at larger interaction; it misses molecular band 
splitting and hybridization with lower 2-particle continuum 
bands and fails to capture the limit of tightly bound molecule 
discussed below. 

At stronger interactions (|.gi_D/a||| ^ hujnm), an accu- 
rate description of the two-body physics requires a number 
of bands of the order 0{huJnm/ Ja,n)- The bound states 
formed below each two-particle continuum at weak interac- 
tions moves downward as the interaction strength increases 
and eventually cross (and in the process, hybridize with) the 
lower two-particle continua (see Figs. [T] (a) and [3] and [13]). 
The specifics of how it hybridizes depends on which band it 
came from [compare the (1,0) curves at A = —3.9 and the 
(2,0) curves at A = —6 in Fig. [T](a), and Figs. [3] (a) and (b)] 
and constitutes the subject of our study. Close to the reso- 
nance, the bound state dispersion relation is strongly modified 
changing the sign J^^-^ [see Fig.[T|(c)]. Figure [ijshows the 
spectrum for the (1,0) and (2,0) lattice resonance deep in the 
tight-binding regime. States inside i?Q q allow the determi- 
nation of the scattering properties of the atoms in the low- 
est band. Also, resonant effects lead to the appearance of 
new bound states at either the edges [Fig. |3](a)] or the cen- 
ter [Fig.[3](b)] of the Brillouin zone (BZ). As shown below, 
the qualitative differences between the (1,0) and (2,0) lattice 
resonances are accurately captured by the if-dependence of 
the atom-dimer coupling. 

Finally, for sufficiently strong interactions (|gi£)/a||| ^ 
hijjnm), the tightly bound dimers which are not in resonance 
with the scattering continuum bands are well described by a 
particle with mass 2m moving in a periodic potential with 
depth 2Vo. Our numerical calculations reproduce this limit- 
ing behavior. 

To gain more physical intuition on the lattice resonance 
phenomena, we adopt a bosonic variant of an effective two- 
channel lattice Hamiltonian Q. For simplicity, we focus on 
the resonant two-body physics in an energy window around 
the (no, no) two-particle continuum. More general (no, mo) 
continuum bands can also be analyzed in a similar manner 
We consider the interaction regime in which the dimer is 
close in energy to the (no, no) atomic continuum band, i.e., 
kd - 2ea,no I < |ea,no+i " ^a,no I whcrc 6^ IS the onsite dimer 
energy. The effective description explicitly introduces a local- 
ized bare dimer state \di) = d\\Q) as the only energetically 
allowed double occupation of a lattice site. The energies of n- 
occupied sites (n > 2) are significantly modified by the strong 
interactions and are expected to fall outside the effective de- 
scription energy window. Thus, we expect that the n-occupied 
sites are energetically suppressed and we impose this assump- 
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tion by treating the atoms and dimers as hard-core objects. 
The exact many-body Hamiltonian can be formally reduced 
to the Hilbert space of empty, singly and doubly occupied 
states [4 | by applying the projector operator Vad = ^arf, i 
where VadA = |0i)(0i| + \ai){ai\ + \di){di\ (here Ui is the 
atomic annihilation operator). The leading interactions terms 
in the tight-binding regime come from nearest neighbor cou- 
plings which involve tunneling of only one atom: i) the cou- 
pling of two atoms to a dimer and ii) the exchange of an atom 
and a dimer on nearest neighbor sites. Keeping only these 
terms, the effective Hamiltonian reduces to 

+9ex ^ d\dja}^ai + ^ gi-j[d\ajai + a\a]di] Vad, 

(3) 

where Ja = Ja,no (Jd) is the atomic (dimer) tunneling, g±i is 
the atom-dimer coupling and g^x is the atom-dimer exchange 
coupling. The Hamiltonian (|3]l can be equivalently under- 
stood as arising from the lattice projection of the two-channel 
model, restricted to nearest neighbor couplings between the 
dominant near resonant closed-channel molecule and nearby 
two-atom continuum. 
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FIG. 2: (Color online) Effective Hamiltonian matrix elements (ME) 
as a function of Vq. 

A crucial aspect of the proposed effective Hamiltonian is 
the dependence of the atom-dimer coupling on the relative 
position {i — j = ±1) between atoms and dimers, g_i — 
(— l)^'*(7+i. This dependence comes from the parity prop- 
erties of dimer Wannier function which are directly related 
to dimer symmetry with respect to the center-of-mass coordi- 
nate. Since the one- and two-body Hamiltonians are invariant 
under a reflection TZ that takes {ri} — > {—ri}, the atom and 
dimer Wannier functions are either symmetric (s-orbital) or 
antisymmetric (p-orbital) under the reflection at the bottom of 
the lattice site in which they are centered. By applying the re- 
flection operation 72. to 5-1-1, and using parity properties of the 
Wannier functions we obtain g_i = {—l)^''g^-i. This struc- 
ture is crucial to capture the hybridization of different excited 



dimers with the two-atom continua (see Fig. [3]). 

We next turn to the quantitative determination of i?e// Pa- 
rameters. This requires a full lattice solution, not just single 
site approximation |4j. Associated to a dimer state, there is a 
dimer Wannier function Wm {ri , r2 ) that can be used to obtain 
the parameters of the effective Hamiltonian. The dimer Wan- 
nier function can be determined by first solving the two-body 
problem in a reduced Hamiltonian H' = V^HV and then us- 
ing the dimer Bloch functions 4>m,K{'''i,r2) to construct the 
dimer Wannier function. Here V projects the Hilbert space 
outside the (no,rio) scattering continuum band, ensuring the 
orthogonality between the dimer and scattering states of the 
no band [V = 1 - V2a with V2a = Y^i^j \aiaj){aiajW. A 
similar procedure has been used to extract the "closed chan- 
nel" dimer in confinement induced resonances 1 14|. 

The reduced two-body Hamiltonian describes all dimer 
bands and scattering continuum bands different from 
(no, no). With the dimer Wannier function (center at 
site j) Wmj{ri,r2) = Wmiri - ja,r2 - ja) = 

^2) in hand, we extract the two- 
body parameters of the effective Hamiltonian by requiring 
that the matrix elements of the H^f / match those of the ex- 
act Hamiltonian, i.e., Jd = —{dj\H\d^j^^), = — 

and g±i — {dj\H\a^jaljj.i). Alternatively, Jd and can be 
estimated from the dimer band assuming that the dimer dis- 
persion relation follows the tight binding prediction ed{K) = 
Ed — 2JdCos{Ka). The g^x coupling can be obtained sim- 
ilarly and is expected to be of the order of the Ja, however 
it only plays a role in a systems with more than two atoms. 
Consequently, we leave its quantitative determination to fu- 
ture publications [11] . 

In our numerical implementation, it is not simple to re- 
duce the exact Hamiltonian to H', so we apply an alternative 
method to extract the dimer Bloch functions. First, we solve 
the exact Hamiltonian at an interaction for which the dimer is 
close to the scattering continuum but still weakly coupled to it. 
Then, we Gram-Schmidt orthogonalize the dimer state from 
the uncoupled scattering continuum (described by hard core 
bosons) to obtain approximate description of the bare dimer 
Bloch function and dispersion relation. Finally, we construct 
the Wannier functions and extract the effective Hamiltonian 
parameters. These parameters, illustrated in Fig. [2] converge 
fast with the lattice size and only a few sites are needed to 
reach convergence. 

The bound state and scattering two-body solutions of the 
effective Hamiltonian ([3]),([T]i can be solved analytically for 
an infinite lattice. When solving the two-body problem at 
fixed K, the atom-dimer coupling is proportional to gs(K) — 
2gcos{Ka/2) for P,„ = and gA{K) = 2gism{Ka/2) 
for P„i = 1. Thus, the atom-dimer coupling is maximum 
at K = Q for s-s coupling but at K = vr/a for s-p cou- 
pling. The two-body bound state energies are given hy E = 
J{K){a + l/a) where J{K) — 2 Jcos(fca/2) and a is the so- 
lutionofJ{K)~ediK)a+{\g^{K)\^~J^{K))a^/J{K) = 
with the constraint \a\ < 1. The analytical solutions allow 
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FIG. 3: (Color online) Two-body spectrum at the (1,0) (a) and (2,0) 
(b) lattice resonances. Comparison of the exact and the model Hamil- 
tonian predictions for a lattice of 21 sites with periodic boundary 
conditions. Calculations for Vo = 20-Er at interactions for which 
the (1,0) and (2,0) dimers are in resonance with the (0,0) scattering 
continuum. Shaded region corresponds to the (0,0) scattering con- 
tinuum. Larger and darker (blue online) symbols correspond to the 
exact Hamiltonian predictions while smaller and brighter (orange on- 
line) correspond to the model predictions. The two results are almost 
indistingui shable . 



from zero to two bound states as found in the exact numer- 
ical solution illustrated in Figs. [T] [3] A second bound state 
appear for the K values at which the coupling is maximum, 
i.e, edges of the BZ for ~ 1 [see e.g. Fig. [3] (a)] and 
center of the BZ for = [see e.g. Fig. [3](b)]. If the 
dimer is inside the scattering continuum and the coupling is 
weak, the system supports metastable states whose lifetime 
can be obtained by analytically continuing the bound state 
solutions into the complex plane [see e.g. small K region 
in Fig. |3](a)]. The analytical solutions also reproduce the 
two-body bound states of the single-band Hubbard model US] 
when \di) = a\a\ |0) / y/2 and the coupling is of the symmetric 
parity type (i.e., gs{K)). Under these conditions, g = V2Ja, 
Cd is the onsite interaction energy U, Jd = 0, and the bound 
state energies are E = sign(C/)^t/2 + 16J2 cos2(fca/2). As 
shown in Fig. |2j the atom-dimer coupling is even stronger 
for excited dimers. For Vq = AQEr, « 1.72Ja and 

g(2,o) l.goJa and these values only change by a few per- 
cents in the range Vq ~ 4-40i?r and are in good qualitative 
agreement with the predictions of Ref. [6 |. 

To confirm the validity of the effective Hamiltonian ([TJ, 
(|3]l,we solve the two-body problem for a finite lattice with 
periodic boundary conditions using both H and H^f / in the 
resonant regime. We found excellent agreement between the 
exact and the effective Hamiltonian predictions [see Figure 

i- 

The s-s and s-p orbital symmetry of the coupling can be 
experimentally probed by analyzing the quasimomentum de- 
pendence of the molecules fraction after a magnetic field ramp 
through a lattice induce resonance. Initially dimers are formed 
in excited bands, which can be achieved by populating atoms 



in excited bands jTSl. Then, interactions are tuned through a 
lattice resonance, and finally, the dimer fraction is measured 
as a function of the dimer quasimomentum. At the two-body 
level, the final dimer fraction is well described by a Landau- 
Zener probability exp(— (5lz) with a Landau-Zener parameter 
Slz \gp{K)\'^ /\a\ where a is the speed of the ramp. Thus, 
the final molecule probability will be mainly affected at the 
center (s-s) or the edges (s-p) of the BZ depending on the 
symmetry of the coupling. 

Our predictions on the general structure of the atom-dimer 
coupling are based only on symmetry properties of the atom 
and dimer Wannier functions and can be easily extended to 
multi-component and higher dimensional systems. The struc- 
ture of the lattice induce resonances in higher dimension 
would be determined by the dimer symmetry in each lattice 
axis direction. Thus, we expect a rich variety of physical phe- 
nomena when the system is close to a lattice-induced reso- 
nance, that extends to higher dimensions. 
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